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ABSTRACT 

We study the dynamics of stellar wind from one of the bodies in the binary system, 
where the other body interacts only gravitationally. We focus on following three issues: 
(i) we explore the origin of observed periodic variations in maser intensity; (ii) we 
address the nature of bipolar molecular outflows; and (iii) we show generation of 
baroclinicity in the same model setup. From direct numerical simulations and further 
numerical modelling, we And that the maser intensity along a given line of sight varies 
periodically due to periodic modulation of material density. This modulation period 
is of the order of the binary period. Another feature of this model is that the velocity 
structure of the flow remains unchanged with time in late stages of wind evolution. 
Therefore the location of the masing spot along the chosen sightline stays at the same 
spatial location, thus naturally explaining the observational fact. This also gives an 
appearance of bipolar nature in the standard position-velocity diagram, as has been 
observed in a number of molecular outflows. Remarkably, we also And the generation 
of baroclinicity in the flow around binary system, offering another site where the seed 
magnetic flelds could possibly be generated due to the Biermann battery mechanisms, 
within galaxies. 
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1 INTRODUCTION 

Periodic variabilities in the maser intensities have 
now been observed in a number of sources, with 


periods ranging 

from few days to several years 

(iGoedhart, Gavlard & van der Walt |2004 iGoedhart et al.l 

2004 I 2 OO 9 I: va 

a der Walt. Goedhart & GavlardI l2009l: 

Arava et ahlboiol: 

Szvmczak et al.ll201llj. Such variabilities 


appear to be charactersitic of the methanol masers, which 
trace the early e volutionary stages of massive star formation 
(|Elhngsenl l2QQ6l l . Massive stars form in regions which are 
deeply embedded in the molecular clouds, thus posing many 
observational challenges. Methanol masers provide direct 
access to these regions and therefore prove to be important 
tool to study the star formation mechanisms, while also 
giving valuable informations on variable conditions in the 
environment where these massive stars are born. This 
provides sufhcient motivation to understand the cause of 
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maser variabilities, revealed by the long time monitoring of 
the methanol masers. 


iNorris et al.l (|l998h have found from their high- 
resolution imaging that some of these sources show arc-like 
maser spots with velocity gradients close to those given by 
Keplerian profile, and thus they argue that the methanol 
masers could be possible tracers of circumstellar discs. But 
we also note that there are some sources which do n ot have 
such linear or arc-like structures (I Walsh ^Ugjj. I t has 

also been suggested bv iMinier. Booth Conw^ ( 2002h that 
some of these methanol masers might be associated with 
molecular outflows or expanding HII regions, and thus might 
originate from regions spatially far from those of the stellar 
objects. 


Much of our understanding of the regions where the 
young stellar objects (YSOs) form in the molecular clouds 
is due to the maser emission from such locations. High ve¬ 
locity outflows with velocities ranging from few kms“^ to 
few hundred kms“^ have been observed, and it is gener¬ 
ally agreed that such outflows occur around YSOs, driven 
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due to stron g stel l y winds in the ear ly evolutionary stage s 
(see reviews ISnel 1 (ll98.tl:lLadal (Il985h: IWelch et, alJ (Il98,^ : 
IShu et ^ (ll987l l: lBachiller et alJ (llQQfil ll. These outflows are 
mostly bipolar and are ubiquitious. Sinc e the discovery of 
these bipolar outflows (|Snell et al.lll980h various attempts 
have been made to understand the physical nature of these 
phenomena. It is generally believed that the collimation has 
to happen very near to the star forming region and it seems 
unlikely that it happens either due to anisotropic environ¬ 
ment far away from the ce ntral region or also due to local 
direction of magnetic fields (ICanto et al.|[l9^ : lKonigl|[l982l : 


iTorrelles et ahlflQSSl : iHever et al.l 


Given the widely accepted view that most stars are part 
of multiple star systems, mostly binaries, and the compelling 
evidence from various maser observations that the phenom¬ 
ena, such as, the maser variabilities, or the collimation of 
outflows, should occur locally, we find it important to study 
the dynamics of stellar wind in the binary potential. This 
also leads us to an important question concerning the gen¬ 
eration of seed magnetic field from initially zero magnetic 
field, and we try to study whether such a possibility exists 
in the same model setup. The standard paradigm to explain 
the orgin of cosmic magnetism involves, first, the genera¬ 
tion of seed magnetic fi elds, which is later amplified due to 
the turbulent dyna mos ([Brandenburg Subramanianll200^ : 
ISubramanianllioOsl i . The presence of baroclinicity is known 
to give rise to the vorticity, or the magnetic field in the elec¬ 
trically conducting plasma, s tarting from zero initi al fields 
(see e.g. lSubramanianI (|2008l i: lModestov et al.l (|2014l B. Here, 
our interest is to only explore if the baroclinicity develops in 
our model system, and any furthur studies focussing on the 
vorticity or the magnetic fields will be taken up elsewhere. 

In the present work, by considering stellar wind from 
one of the bodies in the binary system, we focus on follow¬ 
ing thee issues: (i) the periodic variations in maser intensity; 
(ii) bipolar outflows; and (hi) the generation of baroclinic¬ 
ity. In § 2 we describe our model setup and derive some 
general principles in the rotating reference frame. Results 
from direct numerical simulations (DNS) are presented and 
discussed in § 3. Based on our DNS results, we perform fur¬ 
ther simulations to study maser variability and the issue of 
bipolar outflows in § 4. We conclude in § 5. 


2 THE MODEL 

We present a simple model which involves basics of the two- 
body and the three-body problems in classical mechanics. 
We repeat some of the well known things for completeness 
and the details can be foun d in any standard textbook o n 
classic a l mec hanics (see e. g. IValtonen KarttunenI (|2005l h 
iMorinI tosh T 

Consider a binary system consisting of two bodies, S 
and P, which are rotating around their common center of 
mass, O, in a plane. Let XY Z he the inertial (fixed) coor¬ 
dinate frame in which the two bodies lie in the X T—plane 
with angular velocity vector in Z—direction. Let XY Z he 
the rotating (comoving) coordinate frame which rotates with 
an angular velocity same as that of the two bodies in the bi¬ 
nary system and therefore both the bodies appear to be at 
rest in this frame. The origins of both the coordinate frames 
coincide and are taken to be at the center of mass, O, of 


Y 



Figure 1. The inertial {XY) and rotating (XY) coordinate 
frames are shown for the three-body problem. The center of mass 
of the two bodies labelled as S and P is the origin of both the 
coordinate frames and is denoted by O. 


the binary system. All the assumptions made while study¬ 
ing this problem analytically are given below, but it should 
be noted that while performing numerical simulations for 
this problem, most of the following assumptions are lifted 
and therefore we simulate a more realistic case. 

(i) One of the bodies, S, has a spherically symmetric wind 
very near to its upper atmosphere, whereas the other body, 
P, is interacting only gravitationally. 

(ii) S and P, which are called the primaries, move in cir¬ 
cular orbits around their common center of mass, O, in a 
plane. 

(hi) The mass of the wind is assumed to be negligible as 
compared to the total mass of S and P. Thus, it is essentially 
a restrieted eireular three-body problem, with only difference 
that the wind (massless), which mimics the third body, is 
modelled as a eontinuum. 

(iv) Molecular viscosity of the wind is assumed to be neg¬ 
ligible. 

(v) The wind flow is assumed to be in the steady state. 

(vi) We assume that the fluid behaves as a perfect gas 
and the flow is isentropie. We use perfect gas equation of 
state for our analysis. 

The units of various quantities are chosen such that the 
properties of the system depend only on a single parame¬ 
ter. Let the total mass {M) of the primaries {S and P) be 
the unit of mass; the distance between them (V) be the unit 
of distance; and the unit of time be chosen in such a way 
that the angular speed of the primaries, denoted by Q, be 
unity. It is known that the two bodies circling around each 
other satisfy the following relation: 

p3 ’ 

where G is Newton’s gravitational constant. We find it useful 
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M/Mq 

Figure 2. Relation between the total mass (Ji) and the time- 
period (ST) of the binary. Different lines in this log-log plot corre¬ 
spond to different values of the binary separation {T>). 



- 2-10 1 2 


Figure 3. Isocontours of the effective potential given by 
Eq (O, for ^ = 0.2. Locations of Lagrangian and the temperature 
points are shown by Li — L 5 and Ti — T 5 , respectively. 


to express Eq. © in the following form: 


gr 

IDay 


2 

= 1.33376 X 10® 




( 2 ) 


where & = 27r/Q is the time-period of the binary and Mq 
represents the solar mass. 

In dimensionless units, equation o implies G = 1 due 
to the choices made above for the units of mass, length and 
time; also, the mean anomaly equals time (r = r is the time 
as seen in both the reference frames). Figure 1 shows the 
binary system at an arbitrary time r as seen in both the 
coordinate systems, in which A is an arbitrary point, and 
I* S'A, rpA and roA are the radius vectors of point A with 
respect to S, P and O respectively (the axis Z = Z is not 
shown explicitly). Let ^ be the mass of P, thus mass of S 
is (1 — (f). In the rotating reference frame, with positive X 
in the direction of the body P, the coordinates of P and S 
will be (1 — ^,0) and (—^,0) respectively. Let (X, E, 0) be 
the coordinate of the arbitrary point A, assumed to lie in 
the plane of the binary. Thus magnitudes of vectors rsA and 
rpA are given as: 


rsA = |rsA| = V{X + 0 ^ + Y^ 


rpA = \rpA\ = ViX - (1 - Oy + Y^ (3) 

The gravitational potential in the comoving frame at the 
point A may be written as, 

$ = -iLiil _ (4) 

rsA rpA 


2.1 Bernoulli’s principle in the rotating 
coordinate system 

If u(X,t) be the fluid velocity of the wind in the rotating 
frame then we may write the Euler equations in rotating 
frame for steady flow, with p and p as the fluid pressure and 
density, respectively, as, 

(i;.v)v = -2nxv (5) 


where (X, r) = (X, Y, Z, r) and ft (= ez, which is the unit 
vector along Z = Z) is the angular velocity of the comoving 
frame relative to the inertial frame. Note that the equa¬ 
tion © is written in dimensionless form. Using vector iden¬ 
tities, we can write equation ([5]) in the following form: 

VP + (cp + 2fi) X u = 0 (6) 

where, 

B = Qv" + / y + -^‘eff) ; w = V X V (7) 
The effective potential ($eff) in equation ([7|) is written as, 

= $-i|nxV|" = -ILA)-i-(8) 

2 TsA rpA 2 

On taking the dot product of equation with V, we obtain 
(^;.V)P = 0 (9) 

Therefore the quantity, P, is a constant along a particular 
streamline for steady flows, although it could be a different 
constant for different streamlines. Noting the fact that the 
particle paths and streamlines are the same for steady flows, 
we can see that P remains the same for a particular fluid 
element as it moves along a particular streamline. 


2.2 Nature of the flow 

The Lagrangian time derivative is the rate at which some 
quantity of interest changes as we follow any particular fluid 
element, which is denoted as, 

A.| + ,..V) (10) 

As the first term on the right hand side of equation (nni 
(which is the Eulerian time derivative) does not contribute 
for steady flows, we can write from equations 0 and (IIQ}, 

^ = {vV)B = 0 (11) 
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4 Singh & Deshpande 

Equation m implies that, for any streamline, which is same 
as the trajectory of a particular fluid element, we may write 

B = ^RT $eff = constant = C (12) 

where we have used adiabatic equation of state and note 
that the term J dpip appearing in equation © may be re¬ 
placed by specific enthalpy (le) for isentropic evolution of 
fluid element. We know, 

“ = (7^) <“) 

where 7 is ratio of specific heats at constant pressure and 
constant volume, R is the gas constant and T is the tem¬ 
perature. As the terms u^/2 and { 7/(7 — 1)} RT in equa¬ 
tion (O cannot be negative, we infer from equation cu 
that the motion of a fluid element, and hence the corre¬ 
sponding streamline, is restricted to the region where 

^eff < C (14) 

2.3 Isocontours of effective potential and 
temperature: zero velocity curves 

It is known that the isocontours of the effective potential 
exhibit five Lagrangian points, at which the test-body may 
remain at rest in the comoving frame. These points are ex¬ 
trema or saddle points of 4>eff. It is remarkable to note that 
the shape of the isoeontours of the temperature distribution 
for binary system under eonsideration is identieal to the iso¬ 
eontours of the effeetive potential^ as may be seen from the 
following discussion. 

Let us focus on the zero-velocity-curves, the topology of 
which depends on the energy of a particular fluid element. 
Let $eff = C'l, where Ci is constant. Then, we can see from 
equation for zero-velocity-curves, that 

^^ ~ ~ C — Cl — constant 

^ T — constant, where, 4>eff — Ci ^ (15) 

as { 7/(7 — 1 )}^ is a constant for a particular value of 7 . 
Thus we conclude that the shapes of the isocontours of the 
effective potential and the temperature are same, whereas 
the values might be different. Hence we will have five special 
temperature points^ Ti — T5, at the same locations where we 
have five special Lagrangian points, Li — L5. We may refer 
to these five temperature points as Lagrangian-equivalent- 
temperature points. 


3 DIRECT NUMERICAL SIMULATIONS 

So far, we have studied the behavior of few scalar variables 
along streamlines. We now wish to know the trajectories of 
fluid elements, i.e., we would like to study how the spheri¬ 
cally symmetric wind from the body S flows in the presence 
of another gravitating body P in a binary system. For such 
investigations, we use PLUTO codeQ, which is a Godunov- 
type modular code intended primarily for computational as¬ 
trophysics and high mach number flows in multiple spatial 

^ See http://plutocode.ph.unito.it/ 


dimensions. The code allows us to study two-dimensional 
problems in a frame rotating with constant angular veloc¬ 
ity pointing along Z-direction, by suitably adding nonin- 
ertial (the coriolis and the centrifugal) forces in the mo- 
mentum equat i on. T he details of the code may be found in 
iMignone et al.l (|2007l ) (and references therein). Our strategy 
to study the dynamics of the wind in a binary system using 
PLUTO code may be expressed as follows: 

(i) We adapt the code in the comoving frame of the binary 
in which the two bodies, S and P, appear to be at rest, 
by adding the necessary body-forces, namely, coriolis and 
centrifugal, to the equation of motion. 

(ii) We use the hydrodynamic module of the code and 
solve the equations in two-dimensional r-(f plane, where X — 
r cos 4> and Y — r sin 0, and the angular velocity points in 
the Z-direction. Linearized Roe Riemann solver has been 
used for flux computation and we used outflow boundary 
conditions at the outer edge of the domain. 

(iii) As we wish to study how the wind from one of the 
bodies {S) in binary system flows due to purely gravitational 
effects, we do not consider the effect of forces that might ac¬ 
celerate the wind radially outwards from S when it leaves 
the surface of the body S (e. g. radiation pressure on the 
wind due to S may lead to radially outward acceleration 
of the wind). One may wish to imagine a radially-outward- 
acceleration-zone around body S, beyond which, the dynam¬ 
ics of the wind is solely governed by the gravitational effects 
due to binary system and the pressure gradients. The radius 
of such acceleration zone was chosen to be equal to 20 % of 
the binary separation. 

(iv) Thus we consider the spherically symmetric wind 
from S and study its dynamics in the plane of the binary. 
Equation of motion for the fluid particles is given by: 

^ + (vV) V = - V$eff - 20 X v (16) 

or p 

where symbols have ususal meanings described earlier. As 
our aim is to illustrate the physical mechanism, we are not 
much interested in absolute values of various physical quan¬ 
tities, but their relative changes as one moves about in space 
at any given time, would indeed be useful for our further 
modelling. 

3.1 Spiral nature of outflowing gas 

We present results from two different simulations with ^ = 
0.1 (model A) and 0.3 (model B), and show the evolution of 
density/pressure maps with time as seen in the corotating 
frame, which is the rest frame of the binary system. The 
domain in r- and (f>- extends from 0 to 4 and 0 to 27r, respec¬ 
tively, with V being the distance unit. The grid resolutions 
were chosen to be 300 x 300 for model A and 400 x 400 for 
model B, in r X 0 domain. We reduced the total mass (M) 
of the binary by factor 4 in simulation with ^ = 0.3 as com¬ 
pared to the one with ^ = 0.1. As noted before, Ji is the 
unit of mass in this work. 

In Figs (HI) and ([5|) we show, respectively, the snap¬ 
shots of the density and pressure, both in logarithmic units 
(log^oP model with ^ = 0.1. Time in¬ 

creases from panel ADI (API) to AD 6 (AP 6 ) and we see the 
development of the spiral pattern of the outflowing mate¬ 
rial. The evolution reaches a steady state in later stages and 
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Figure 4. = 0.1) Snap-shots of logarithmic density maps, as seen from the corotating frame; colours indicate log^^Q p. Time increases 

from upper-left to lower-right panel. Two concentric circles are centered about the common center of mass of the binary, denoted by the 
asterisk (‘*’). The dashed line joins the centers of the two components of the binary, where the star, S (with the wind), lies at the left 
end, whereas the other body, P, lies at the right end. Arrows indicate the velocity vector field. 


© 2012 RAS, MNRAS 000, [THE] 


























6 Singh & Deshpande 




2 


^ 0 


-2 



-L22 

-L77 

-2.31 


- 2.86 


-3.40 


-3.95 


-4.49 


Figure 5. Similar to Fig. but shows the snap shots of the pressure maps on logarithmic scale (log^Q P). 
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Figure 6 . Same as Fig. (|4|) but for ^ = 0.3. 
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Figure 7. Same as Fig. (| 5 | but for ^ = 0 . 3 . 
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thus we find that the structure of the velocity field stops to 
evolve considerably. We note that the symmetry axis of the 
spiral pattern is shifted from the rotation axis (which passes 
through the center of mass of the binary) by distance The 
whole pattern of density/pressure as shown in Fig dH)/® 
rotates with binary period for an inertial observer, and thus 
we find periodic variation in the mass density along a given 
sight line. To illustrate this, we have drawn two concentric 
circles, centered around the common center of mass of the 
binary, showing the density/pressure variations as one goes 
around the circles. The period of such density modulations 
will be of the order of binary period. 

Similarly we show snapshots of density and pressure for 
the model with ^ = 0.3 in Figs © and 0 , respectively. 
As before we find that the stellar wind spirals outward and 
reaches a steady state. The structure of the velocity field 
stays nearly the same in late stages, which would be cru¬ 
cial to explain the observational fact that the maser spot 
does not move in the sky. Although our aim is to propose 
a physical model which can explain some key qualitative 
features seen in the observations related to maser intensity 
variations, we find that a closer look of Figs ©-© reveals 
that the density at some spatial location along a particular 
sight line might vary easily by a factor of about 10 between 
its maximum and minimum values. 

To make contact with physical units, we refer the reader 
to Fig. m where we show the relation between the total 
mass (M) and the time-period ( 21 ) of the binary, corre¬ 
sponding to different values of the binary separation (V); see 
also Eq. ©• The mass, period and distance are expressed 
in terms of solar mass, day and astronomical unit (AU), re¬ 
spectively. If we let the distance unit in our simulations to 
be equal to 1 AU, the periods corresponding to total masses 
equal to 1 M© and 50 M© will be about 1 year and 50 days, 
respectively. 


3.2 Baroclinicity 


We now turn to a completely different phenomena which 
will have implications for the generation of the vorticity, 
given by Cc? = Vxv. Given that the flow around binaries 
in star-forming regions are expected to be highly ionized, 
this is usually studied under standard magnetohydrodynam¬ 
ics (MHD). Although the study of magnetic fields and the 
plasmas is beyond the scope of the present paper, we, how¬ 
ever, wish to remark that the magnetic field, B, could also 
evolve due to the presence of baroclinicity, in the same model 
setup. The mathematical form of the evolution equations of 
the vorticity and the magnetic field are same; these may be 
expressed as. 


dB 

dt 


Vx(vXB) + 


'VPx'Vpp 

pI 


+ ??v"b 


doj 

dt 


VX (vxu?) 


VPx Vp 


P 


2 




(17) 

(18) 


where 77 , u and pp are the magnetic diffusivity, the kinematic 
viscosity and the plasma density, respectively. The second 
term on the right hand side of Eq (HU (or HU is known 
as the baroclinic term, which can lead to the generation 
of B (or u?) even if the magnetic field (or vorticity) was 
strictly zero to begin with. The baroclinicity characterizing 


the misalignment between the isocontours of the pressure 
and the density may be expressed as. 


0 = 


VPx Vp 

p2 


(19) 


In the same model being discussed here, the baroclinic 
vector field /3 is zero everywhere initially. It is instructive 
to see if such a field could be generated at later stages of 
evolving stellar wind in the Roche potential. To study this 
we considered both the models A and B, discussed above, 
for which we have two dimensional data (in the plane of the 
binary) of the density and the pressure fields at different 
instances of times; few snapshots of density and pressure 
are shown in Eigs ©-©• Erom this data, we compute /3 
using Eq (1191) which points along ±ez by definition, as both 
the vector fields, VP and Vp, lie in the plane of the binary 
in our 2D setup. We find that (3 changes sign in the binary 
plane, but its magnitude is quite large, especially along the 
spirals. In Eig. (| 8 ]) we show the snapshots of the quantity 
log^o l/^l models A (top panels) and B (bottom panels), 
in the late stages of the simulations. 

Thus in the same setup, we show that the baroclinic¬ 
ity develops significantly, which can source the vorticity field 
and also the magnetic field if the fluid is electrically conduct¬ 
ing, according to Eqs. (HB and m- Similar conditions are 
expected to exist in many parts of a galaxy, and therefore 
this provides yet another possibility by which the seed mag¬ 
netic fields could arise within the galaxy. The spatial scales 
of variations of /3 are smaller compared to the binary orbit 
and therefore this could source the seed magnetic (or vortic¬ 
ity) fields at small scales. Detailed investigations focussing 
on the magnetic fields, giving quantitative estimates, are 
needed, and will be studied in a future work. 


4 MASER VARIABILITY (LIGHT-CURVE) 

Now we focus on the possible distribution of maser sources 
in the medium around the binary system, as seen by the ob¬ 
server in the defined geometry. The apparent (line-of-sight) 
velocity structure and the column density of the molecular 
matter (for relevant species) are amongst the key ingredients 
governing the formation of astronomical masers. Hydrody¬ 
namic simulations presented in § 3 reveal the flow structure 
and distribution of column density in our model setup. We 
make use of these essential informations and perform further 
numerical simulations to study the masing action in such 
environment. Eor simplicity, we have used a single-arm spi¬ 
ral pattern in these simulations, instead of a two-arm spiral 
structure as suggested from PLUTO simulations presented 
in the previous section. We note that the patterns of field 
variables (mass density, pressure and baroclinicity), shown 
in Eigs (ll])-®, rotate with the angular speed (U) of the 
binary, for an observer in a fixed inertial frame. In this sec¬ 
tion we present our results as would be seen by an inertial 
observer. 

Physical mechanisms governing the laws of maser for¬ 
mation in cosmic se ttings are complex and have been studied 
extens i vely (see e.g.lGoldreich KeelevI (|l972l i ; IStrernitskiil 
(ll974h : lElit,zn3(ll992ll l. A widely accepted i nodel for class II 
methanol masers is bv ISobolev Deguchil (|l994li . Here our 
focus is on the conditions which can potentially cause peri¬ 
odic variations in maser intensities, while also being able to 
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Figure 8. Snapshots of the quantity, log^Q |/3| (shown by colours; |/3| = |/3z| in this setup), for models A (top panels; ^ = 0.1) and B 
(bottom panels; ^ = 0.3). Left and right panels correspond to panels 4 and 5 of Figs. (|4l)-([7|). 


address the observational fact that the sky locations of the 
maser spots stay fixed with time. We do not intend to study 
the formation mechanisms of maser and refer the reader to 
earlier works, some of them quoted above, for further related 
studies. For our purposes, it is sufficient to note that the in¬ 
tensity of maser, Z oc exp(—r), where r is optical depth 
(negative for amplification), which varies along a particular 
sight line, and can be expressed in an average sense as. 


where E is the mass column density of the medium and Ai’y 
is the velocity spread along the chosen line of sight. Thus 
the sightlines with vanishing velocity gradients are ideal for 
masing action. If the structure of the velocity field in the 
sky does not change, then the masing spot is expected to 
stay at the same sky location. 

Let IZ denote the location of the masing spot in the 
sky, seen from the fixed coordinate frame XYZ. If the 
density of the medium at IZ varies sinusoidally in time as, 
pilZ, t) ~ po sin {cut + (/?)+ yOf, then, to lowest order in r, the 
maser intensity will also show sinusoidal temporal modula¬ 


tion with the same period (= 2tt as may be seen from 
Eq (1^ . Thus any temporal modulation in the column den¬ 
sity will result in maser intensity variation. The quantities 
Pq^ (f and pf denote, respectively, the amplitude, angular 
frequency of modulation, the phase and the floor density. 
Given that the density maps shown in § 3 rotate with an¬ 
gular speed (Q) of the binary, for an inertial observer along 
a chosen line of sight, the observed density at IZ will show 
modulations with the angular frequency (cl;) which is com¬ 
parable to Q, i.e., uj ~ 0(0). As we have modelled a single¬ 
arm spiral structure, for simplicity, in this section (instead 
of two-arm spiral structure suggested from PLUTO simu¬ 
lations), we note that the orbital periods will be overesti¬ 
mated by factor two. We emphasize again that our aim here 
is to explore possible physical conditions that can lead to 
observed periodic modulations in maser intensities. One of 
the conclusions from our study may be stated as: the binary 
period determines the period of maser intensity variations. 

Binary orbits might be inclined at arbitrary angles with 
respect to sight lines. The inclination angle (za) may be de¬ 
fined as the angle between normal to the binary plane and 
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Figure 9. Maps of (i) column density (top left); (ii) distance-velocity (top right; colours: maser intensity); (iii) component of velocity 
along the line of sight (bottom left); and (iv) the gradient of velocity along sightline (bottom right), in a fixed inertial frame, shown for 
a binary system with angle of inclination = 60°. Colors indicate the corresponding quantities on arbitrary scale. 


outwardly pointed line of sight. We wish to illustrate the 
emission patterns when the binary lies somewhere between 
the edge-on (za = 90°) and the face-on (za = 0°) configura¬ 
tions. Here we consider a binary with za = 60°. We identify 
maser emitting regions in our simulations and study the tem¬ 
poral evolutions of the maser intensities from many different 
sources. 

Figure Q shows maps of various physical quantities of 
interest at some time and we briefly discuss these below: 

(i) The column density shows single-arm spiral structure 
where the density decreases as we move away from the cen¬ 
tral regions; top left panel. 

(ii) Remarkably, the distance-velocity diagram in top 
right panel shows clearly the bipolar structure, as has been 
widely observed. This is a natural outcome of our studies 


which suggests that the observed bipolar nature of molecular 
outflows could also arise due to presence of binary systems. 

(iii) Bottom left panel shows the map of the component 
of velocity along the line of sight. It is useful to note that the 
structure of the velocity field remains the same with time in 
this setup and therefore the potential maser spots do not 
move in the sky, thus explaining this observational fact. 

(iv) The map of velocity gradient along the line of sight 
is shown in bottom right panel. Low gradients in velocity 
determine the potential sites where the masing action can 
take place, as may be seen from Eq (|^. 

This leads us to construct and monitor maser inten¬ 
sity maps as functions of time. In Fig. (US we show a typ¬ 
ical map of maser intensity at some time (top left panel), 
where we mark the locations of five different sources (masing 
spots). We intend to show the variations of maser intensities 
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Figure 10. Maser intensity map around a binary system with an inclination angle za = 60° shown in top left panel, where locations of 
five masing spots are marked; seen from a fixed inertial frame. Lightcurves of all these five sources are shown in the rest of the panels. 


as functions of time, or equivalently, as functions of orbital 
cycle of the binary. Such plots of intensity variation, com¬ 
monly called light-curves, corresponding to all five sources 
are shown in other panels of Fig. (HOI). Some noteworthy 
properties are as follows: 

(a) We first notice that the maser intensities show peri¬ 
odic temporal modulations. These modulations are caused by 
the periodic variation of column density, as discussed above, 
and the period is determined by the binary period. 

(b) Multiple masing spots could form around a single bi¬ 
nary system, and ah these sources have same periods of mod¬ 
ulation, as they arise due to density variations, the period 
of which, in turn, is determined by the binary system. How¬ 
ever, the amplitudes of intensity variations are different for 
different sources. 

(c) Thus we argue that the observed periods of such in¬ 
tensity variations trace the characteristics of the binary, i.e., 
the binary period. This, together with some other indepen¬ 
dent estimates of either the mass or separation of the binary, 
could be very useful in determining the nature of binary sys¬ 
tems, embedded in molecular clouds. 

(d) Both, the arc-like and somewhat symmetric spots, are 
possible; compare, e.g., spots ‘2’ and ‘3’ in Fig. (ITOl) . This 
might depend on the local physical conditions and the dis¬ 


tances of the spots from central regions. Linear structures 
could indeed be reminiscent of K eplerian-like velocit y gra¬ 
dient, as has been conjectured bv iNorris et al.l (|l998h . 


5 CONCLUSIONS 

With an aim to provide physical mechanisms that could po¬ 
tentially cause observed periodic variations in the maser in¬ 
tensities, we studied the dynamics of stellar wind from one 
of the bodies in the binary system. We find that the intensity 
variations are due to the periodic variations in the local ma¬ 
terial density, where the period is on the order of the binary 
period. The underlying velocity structure is Keplerian-like 
and it remains frozen-in-time in the later stages. We note 
that this non-evolving velocity structure is important for the 
masing spots to stay at the same sky location, as has been 
revealed by the observations, and it is the maser intensity 
which varies periodically. The wind appears to spiral out¬ 
wards before settling into the faraway Keplerian orbits, and 
naturally appears bipolar in the standard position-velocity 
diagram. 

The spherical symmetry of the stellar wind is broken 
due to the influence of the other gravitating body. This sym- 
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metry breaking very near to the star in the binary system 
holds the key to establishing density inhomogenieties also 
far away from the binary location where the conditions for 
masing action are favourable. We first studied this probelm 
in the rotating frame in which the binary components are at 
rest. To an inertial observer, such density inhomogeneities 
in the rotating frame lead to a periodic density variations 
along a chosen line of sight. Sightlines with nearly vanish¬ 
ing velocity gradients along those directions being ideal for 
observing masers show sinusoidal modulation in the mate¬ 
rial density, thus causing the maser intensity also to show 
sinusoidal variation. The structure of the velocity held be¬ 
ing nearly frozen, which is close to Keplerian, explains why 
these periodically varying maser spots do not move. This 
mechanism also naturally gives a bipolar appearance in the 
standard posit ion-velocity diagram. We note that the bipo¬ 
lar outhows are ubiquitously observed in nature and are pre¬ 
sumably thought to be associated with star forming regions 
in molecular clouds. A binary origin may be responsible for 
some of these bipolar outhows. 

In the same setup, we also hnd that the baroclinicity 
develops and it is mostly concentrated in the spiral form; 
see Fig. This offers yet another scenario, which can lead 
to the generation of seed vorticity and magnetic helds within 
the galaxies, and might further affect the evolution of these 
helds. More detailed investigations focussing primarily on 
the magnetic helds will be presented in a future work. 
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